Astronomy &: Astrophysics manuscript no. 0100 


February 2, 2008 


(DOI: will be inserted by hand later) 





Collisions between carbonaceous grains in the interstellar 

medium 

R. Papoular^ 



cn 
o 
o 

> 

o 



> 

oo 



m 

O ' 

Oh 
I 

o 



Service d'Astrophysique and Service de Chimie Moleculaire, CEA Saclay, 91191 Gif-s-Yvette, France 
e-mail: papoular@wanadoo.fr 



Abstract. An abstract should be given 



Semi-empirical molecular dynamics is used to simulate 
hydrocarbon grain sputtering and collisions which are ex- 
tremely difficult to study experimentally. This microscopic 
and dynamic approach is particularly suited to high ve- 
locity impacts, where target destruction occurs far from 
equilibrium. A wide variety of processes are encountered, 
depending on grain size and velocity: vaporization, frag- 
mentation, atomic implantation, sticking, elastic recoil, 
atomic chemisorption, H abstraction and H2 formation, 
etc. The impact-velocity threshold is about 10 km/s but 
complete grain destruction requires much higher velocities 
and nearly equal grain sizes. The main outcome of strong 
collisions is vaporization, i.e. formation of small molecules 
rather than solid fragments. Most of the impact energy is 
carried away by these molecules in the form of kinetic en- 
ergy, with the consequence that the destruction efficiency 
of collisions is considerably reduced. 

Key words. ISM: dust, ISM:kincmatics and dynamics, 
molecular processes 



1. Introduction 

Strong collisions between IS (InterStellar) dust grains 
are commonly expected to alter significantly, if not out- 
right determine, the grain size distribution, mainly by de- 
stroying big grains thereby creating small ones, but per- 
haps also by creating large grains through coagulation 
of smaller ones (see Wickramasinghe 119(571 Bierman and 
Harwitt 198fl Spitzer ri978|l . If grain destruction is pre- 
dominant in the ISM (IS Medium), then, grains larger 
than 1000 A are expected to be almost depleted; since 
they are observed to be quite abundant, the far-reaching 
conclusion may also be drawn, that grains must also be 
able to grow even in such a tenuous environment (see 
Boulanger and Cox 1998: Miville-Deschenes et al. 
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Also, if collisions induce vaporization, new molecules can 
be created, thus contributing to the rich population of ob- 
served molecules in the ISM. Finally, even if shattering 
does not occur upon collision, intense heating may cause 
structural changes in the grains. These expectations are 
significant enough to have triggered a lasting interest in 
grain collisions for at least half a century (see Salpeter 
[Wm Barlow and Sil kl977. CowielMHI Spitzer [llT78i and 
Draine and Salpeter Il979fl and produced a large body of 
literature. Since grains of relevant size are very difficult 
to produce, handle and accelerate in the laboratory, it is 
not surprising that this literature is essentially theoretical 
in nature, although, of course, it does rely upon whatever 
experimental results are available. 

Earlier treatments considered vaporization to be the 
dominant result of grain collisions, and to be complete 
beyond the relative velocity at which the projectile ki- 
netic energy was equal to the total dissociation energy. 
Assuming the collision is totally inelastic and all the ki- 
netic energy is available for dissociation, this velocity is 
independent of grain size and estimated at 3 km/s for 
icy dielectrics up to 8 km/s for iron particles (Spitzer 
[H^ . More recent papers on the subject have relied upon 
condensed-state shock theory in order to predict effects of 
grain collisions (see Borkowski and Dwek I1995L Tielens 
et al. 119941 Seab il990). This approach enhances shat- 
tering and fragmentation over vaporization, with a con- 
sequent lowering of grain destruction threshold down to 
relative grain velocities of about 1 km/s, and complete 
target disruption at 100 km/s. Admittedly, however, such 
treatments apply best to macroscopic bodies and are not 
easily extended to the grains of interest here. Moreover, 
little data is available on shocks in materials similar to 
those that are thought to be the constituents of IS grains. 

The theoretical and experimental obstacles encoun- 
tered in the field of small grain collisions, even outside 
the realm of astronomy, and the development of mate- 
rial and chemical computation software have prompted a 
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flurry of numerical simulations which shed a striking new 
light on the subject (see Cleveland and Landeman 1992] 
Schulte 1995 ; Aoki et al. 1997). Of course, because of lim- 
itations on available computer memory, this approach is 
best suited to relatively small systems, but these are pre- 
cisely those we are interested in. Moreover, the composi- 
tion and structure of the interacting systems may be tai- 
lored to suit the problem at hand. Also, various relevant 
microscopic phenomena can be followed quantitatively in 
time and space with exquisite detail: atomic adsorption, 
abstraction and implantation, molecule formation, grain 
cratering, shattering or vaporization, etc (see Davenas 
I1993|l . Again, by contrast with shock theory, which sub- 
sumes thermal quasi-equilibrium, the time scales of com- 
putations can be short enough (10~^^s) to accomodate the 
out-of-equilibrium phenomena that are expected to occur 
upon high- velocity impacts of small systems. 

These potentialities are taken advantage of to revisit 
grain collisions in the present paper. The scope here is 
limited to carbonaceous material similar to the presently 
envisioned constituent of IS grains: partly aromatic (6- 
and 5-C rings), partly aliphatic (C-H chains) and more or 
less hydrogenated. Their sizes range from 1 to 60 atoms 
and their velocities from 1 to 300 km/s. The penetration 
range of single C and H atoms in this material is first 
determined as a function of impact velocity. This already 
gives some insight into the mechanism and extent of ma- 
terial destruction upon collision. More quantitatively, a 
cross-section for bond dissociation, and the corresponding 
threshold, are determined as a function of impact velocity. 

Impacts of clusters of the same model material are 
then analyzed, which provides the average dissociation en- 
ergy of the material, as well as rudimentary statistics of 
product sizes. These results are finally used to estimate the 
ranges of grain size and velocity in which strong size mod- 
ifications are to be expected. Generally speaking, these 
estimates are less pessimistic than is usually the case for 
the survival of C-rich grains in the ISM. However, no dis- 
cussion is offered here of the detailed effects of collisions 
on grain size distributions. 

2. Designing the collision experiment 

A description of the software package used for simulation 
can be found in Papoular ((2002|l . The particular code used 
here is AMI (Austin Model 1), as proposed by Dewar et 
al. I|1985|l . This semi-empirical method combines a rig- 
orous quantum-mechanical formalism with empirical pa- 
rameters obtained from comparison with experimental re- 
sults. It computes approximate solutions of Schroedinger's 
equation, and uses experimental data only when the Q.M. 
calculations are too difficult or too lengthy. This makes 
it more accurate than poor ab initio methods, and faster 
than any of the latter. 

First, a molecule of the desired composition, structural 
type and size is built on the computer screen by picking 
each atom in turn from a tabic of elements. The code is 
then used to optimize the molecular geometry, i.e. deter- 
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Fig. 1. Set up of a collision simulation. The lone C atom 
(right) is the projectile, with an initial velocity directed 
along the primary axis of the target. Target (left) made 
of amorphous hydrocarbon material: 36 C atoms (larger 
spheres) linked in chains, 5- and 6-membered rings, and 24 
H atoms occupying some of the available carbon valences. 
Length --12 A, height and depth --8 A. 



mine the nearest architecture which minimizes the total 
potential energy of the system. Here, the initial spatial 
distribution of atoms was chosen so that, upon optimiza- 
tion, the structure included short chains as well as 5- and 
6-membered rings of carbon atoms. Hydrogen atoms were 
then added randomly to simulate partial hydrogenation. 
At the end of the process, the distances between atoms are 
indeed those expected for hydrocarbons. The process is re- 
peated to form the second reactant cluster of the collision. 
The cluster sizes adopted here are a compromise between 
computational time and relevance to real IS dust. 

Then, the reactants are set at an initial distance such 
that interaction between them is negligible: more than 3 
A. Finally, one of them, the projectile, is given, as a whole, 
the desired velocity towards the other, the target; the lat- 
ter is given velocity and the calculation is started. In 
the example of fig. 1, the target is a 60-atom cluster and 
the projectile, on its right side, is comprised of a single C 
atom. 

The molecular dynamics relies on the Born- 
Oppenheimer approximation to determine the movements 
of atoms under nuclear and electronic forces due to their 
environment; the calculation is made by steps of 10"'* ps 
and 10"^ ps respectively for high and low relative velocies. 
At every step, all the system parameters are memorized 
as snapshots so that, after completion of the run, a movie 
of the collision can be viewed on the screen, and the tra- 
jectory of any atom followed all along. This makes for a 
better understanding of the details of destruction mech- 
anisms and outcomes. Note that the dynamics of bond 
dissociations can only be simulated by using Unrestricted 
Hartree-Fock (UHF) wave functions in the Q.M. part of 
the calculation (see Szabo and Osthmd il989i|l . 
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3. Atomic carbon projectile 

A projectile made of a single C atom, as in fig.l, is a 
good start since it is easier to follow its trajectory and 
energy after it hits the target. Also, its behaviour already 
encapsulates the main features of grain collisions for, at 
the high velocities of interest here, the translational energy 
of a projectile atom is much larger than its binding energy, 
so that, even if it is part of a cluster, it behaves nearly 
independantly of its neighbours. 

At the beginning of each numerical experiment in this 
series, the target will be set at rest and the projectile, a 
few angstroms away, will be given a velocity Vi towards the 
target. Since gaseous flows in the ISM usually enforce rela- 
tive velocities between grains, Vi will be taken as the main 
control parameter rather than the kinetic energy. Vi will 
be expressed in units of A/ps (0.1 km/s). It is found that, 
for Vi <100 A/ps, the C atom most often simply sticks 
onto the nearest "acceptor" site in the target. This is a 
result of the strong reactivity of atomic carbon, which has 
four dangling (free) bonds, as well as of the reactivity of 
the target, whose carbon bonds are not all saturated with 
hydrogen. In such a sticking event, the target structure 
is hardly modified, except, perhaps, for the displacement 
of one of its hydrogen atoms. Incidentally, this sticking 
process is also an elementary step in the process of grain 
growth. Not being interested, at this time, in events of this 
type, we now focus on the velocity range beyond 100 A/ps 
(10 km/s). Note that at this velocity a C atom carries a 
kinetic energy of 6.24 eV, somewhat higher than the bond 
energy in hydrocarbons. By way of comparison, Salpeter 
H1974|l set the threshold velocity for sputtering by heavy 
atoms at '^30 km/s. 

Figure 2 apphes to the case of Vb=1000 A/ps (100 
km/s). Plotted as a function of time are the kinetic en- 
ergy of the projectile (a), that of the target (c), and the 
variation of the target potential energy (b), throughout 
the "reaction" . The interpretation of these curves is made 
easier if, after the computer run, the accelerated motion 
of the computation is replayed on the screen and snap- 
shots of the system are taken at characteristic times along 
the X-axis of fig. 2. It is found that the projectile hits the 
target at t=0.0025 ps and leaves it at 1=^0.01 ps, after de- 
scribing an irregular trajectory about 7 A in length. Large 
path deviations occur at 0.0039 and 0.0085 ps, when the 
projectile experiences a close encounter with a C atom of 
the target. In such events, the projectile is strongly dece- 
larated but soon recovers most of its velocity because, in 
the present range of velocities, the dominant interaction 
with the target atoms is a "quasi-elastic" recoil of nu- 
clei (see Sigmund ll969|l . Hence the dips in projectile K.E. 
(kinetic energy) . Correspondingly, the target momentarily 
gains P.E. (potential energy): peaks of curves (a) and (b) 
are of opposite signs. However, all the while the projec- 
tile is inside the target, the general trend is a reduction 
of its velocity as it loses energy to bordering target atoms 
in the form of momentum. Hence the initial rise in target 
K.E. (curve c). Part of this gain is readily converted to 
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Fig. 2. Energetics of a collision experiment (configuration 
of fig. 1 with initial projectile velocity=1000 A/ps). a) 
K.E. of the projectile; b) variation of the target P.E.; c) 
(thick line) K.E. of the target. The two pairs of opposite 
narrow peaks correspond to close encounters between the 
projectile and one target C atom. Here, following one such 
collision, the projectile left the target prematurely at 0.01 
ps with velocity 887 A/ps. Exchanges between target K.E. 
and P.E. are still visible afterwards. 



target P.E, leading to a concomitant rise in curve (b). For 
impact parameters such that the energy transferred from 
projectile to target atom is larger than the bond energy, 
part of this energy goes into breaking a bond, and the rest 
into translational and rotational K.E., as well as internal 
energy (vibrational K.E. and P.E. and, possibly, deforma- 
tion or restructuring energy) of the fragments produced. 
Hence the small bumps of opposite signs in curves (b) and 
(c), accompanying the large recoil peaks. 

In the end, only about 30 % of the projectile energy 
was absorbed by the target. This is only due to the acci- 
dental early exit following a close encounter with a heavy 
target atom, and it will be shown below that the stop- 
ping range for this initial velocity is much larger than the 
present target thickness. In the case at hand, internal vi- 
brations are only weakly excited in the fragments, as evi- 
denced by the smoothness of the later parts of curves (b) 
and (c). However the final gain in potential energy is 50 
eV. Most of this was used to dissociate the 6 observed 
fragments. The rest went into structural changes in the 
fragments as compared to the original material. 

Note that the process of fragmentation and deforma- 
tion occured in less than 0.04 ps, i.e. much less than the 
thermalization time (numerically determined for this type 
of material to be longer than 1 ps; see Papoular l2002|l . and 
even shorter than the transit time of a shock wave through 
the target. It is clear, therefore, that thermodynamics can- 
not apply here; indeed, the sputtering phenomenon is usu- 
ally treated at the microscopic level (see Sigmund I19H9I 
and SchiottHnSHIl- 
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Experiments at lower initial projectile velocities show 
that, as Vi decreases, the number of fragments produced 
and the residual projectile velocity when it leaves the tar- 
get both decrease. When Vi <300 A/ps, the projectile 
stops within the target, at a depth which decreases with 
the initial velocity. Also, its trajectory inside the target be- 
comes more and more wrinkled and ends up as a random 
walk until it finally settles in some available "acceptor" 
site (dangling bond or defect site). This, then, is a case of 
atomic implantation. 

In order to quantify this behaviour, we proceed by 
analogy with the usual method, which assumes that the 
loss, dE, of projectile K.E. is proportional to the density 
of target nuclei and to the linear element, dx, of its trajec- 
tory: —dE/E = nadx, where a is a cross-section describ- 
ing the interaction. Here, we replace the denominator E by 
e, a typical bond energy: this ensures that the dimension 
of a is still that of an area and, at the same time, clar- 
ifies the physics, since the left hand side now represents 
the energy loss normalized by a characteristic energy of 
the system. Then, if m is the projectile mass, and V its 
velocity 

- VdV = (ntlm)adx. (1) 

The target mass does not appear here because it is as- 
sumed to be much larger than that of the projectile. The 
results of the numerical simulations, as well as the physics 
of these experiments, are incompatible with a constant 
cross-section. But agreement is good enough if it is as- 
sumed that : (7 = a + hV + cV^. Clearly, the dimensions 
of a, b and c are, respectively, those of an area, an area 
per velocity and an area per velocity square. If not stated 
otherwise, the units of length, area and velocity will be, 
respectively. A, A^ and A/ps. Then, the solution of eq. 1 
is of the form 

xt-Xi = aln )+l3\n{- 2 

m Vi — a V i — p 

where i and f designate initial and final values, and a, 
(3 and 7 (in the same units as the velocity V) are functions 
of a, b and c. We seek values of a, (3 and 7 so that eq. 2 
is best satisfied by the triads, Vi, V{ and Ax — Xf — Xi, 
corresponding to this series of 7 numerical experiments at 
different V^'s: 150, 200, 300, 400, 500, 750 and 1000 A/ps. 
We find: a=1250 , /3=-130 and 7 = 5.610^2 Deriving a,b 
and c from these values of a, (3 and 7, and substituting in 
the quadratic expression for cr, it is finally found that eq. 
2 is satisfied to a standard deviation of 2 % if 

(7 = 6.6 + 5.110"V- 4.110^^2 (3) 

where a is in A^, V in A/ps, and we adopted n — 10^^ 
cm~^, m=12 atomic units and e=5 eV. 

Note that a is of the order of a few atomic radii 
squared, meaning that, in this range of velocities, the pro- 
jectile hardly interacts with target atoms off its trajectory 
by more than 4 A. Also, from eq. 2, an expression can 
be derived for the range of the projectile as a function of 
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Fig. 3. Line R: range of a carbon atom in the model ma- 
terial of fig.l, as a function of its initial K.E., accord- 
ing to eq. 2 and assuming a material density of 2 g/cm^; 
R(/ig/cm^)=0.02 Aa; (A); squares (E): range measured in 
individual simulations at discrete initial velocities. 

its initial velocity. Indeed, when the projectile's trajectory 
ends within the target, its final velocity can be defined by 
the condition that it become random. In the present ex- 
periments, this is found empirically to be of order 50 A/ps, 
which will be adopted for Vf. The range, Aa; is obtained 
by inserting this value in eq. 2. While Ax is a length, ex- 
perimental studies with different materials show that the 
energy loss is less intimately connected with the length of 
the particle trajectory in the target than with the mass 
of target material swept by the projectile along its path, 
which we will denote by R and is equal to p. Ax, where p 
is the density of the target material. If p~2 g/cm"^ and R 
is in /ig/cm^, then R=0.02 Ax. This awkward dimension 
for R is universally used by experimentalists because it en- 
sures that the corresponding numerical values can be han- 
dled more easily, as illustrated by fig. 3, which represents 
the present results in these units, and can therefore be 
compared directly with published experiments. The nearly 
linear functional form as well as the absolute values of 
this curve are in line with previous sputtering studies (see 
Schiott 119681 for instance) and numerical simulations (see 
Cleveland and Landeman ll992l Aoki et al. ll997|l . The dis- 
tances travelled by the projectile inside the target in the 7 
simulations are also plotted as a scatter graph (E) in fig. 
3. Beyond 400 A/ps (K.E.~100 eV), the range according 
to eq. 2 is longer than the long dimension of the target (12 
A) but, of course, the ordinates of "experimental" points 
cannot reach much higher than this value, so they fall be- 
low the curve. 

Because of the discrete nature of the target material 
structure, Aa; and Vf depend slightly (and randomly) on 
the exact direction of the initial projectile velocity, so that, 
even below 100 eV, there is some scattering of points E 
about curve R. Other illustrations of this scattering are 
provided by fig. 8 to 10. Such scattering is hmited if the 
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grain structure is not too anisotropic and if care is taken to 
avoid head-on atomic collisions and subsequent premature 
exit of the projectile from the finite target (cf fig. 2). 



4. Atomic hydrogen projectiles 

Several simulations were made with initial velocities from 
25 to 3000 A/ps (2.5 to 300 km/s) and various impact 
angles on targets of the same structure as above. In none 
of these cases did fragmentation occur. Below 30 A/ps, 
the projectile simply recoiled from the target. Between 
50 and 100 A/ps, the privileged outcome was abstraction 
of a target H to form a hydrogen molecule; between 100 
and 500 A/ps, chemisorption on the target was preferred. 
Finally, at 1000 A/ps and higher velocities, the projectile 
traversed part or all of the target, leaving only a small 
fraction of its energy in the form of target vibrations. 

The maximum projectile energy in these simulations 
was about 500 eV. At this energy, the sputtering yield of 
on graphite was measured to reach a maximum of only 

Thus, H atoms are not 



0.1 (Vietzke and Philipps ITi! 
expected to contribute significantly to the destruction of 
grains; the physical reason for this is that their mass is so 
low (SalpeterUnnil- 



Fig. 4. Snapshot of a collision between two clusters of 
equal sizes (40 atoms) at t=0: projectile (right) with ini- 
tial velocity 100 km/s towards the target (left); distance 
between centers of mass: 11 A. 



5. Cluster projectiles 

Armed with an understanding of the behaviour of single- 
atom projectiles, we can now turn to projectile sizes of the 
same order as that of the target. Figure 4 shows a series 
of snapshots of target and projectile taken during their 
collision. Each is made of 24 C and 12 H atoms, forming 
a structure similar to that of fig. 1; the initial velocity of 
the projectile (right of fig. 4) is 1000 A/ps towards the 
target (left), which is 11 A away . At t=0.01 ps (fig. 5), 
the projectile is entirely inside the target, with no dra- 
matic change observed in either structure. At this stage, 
momentum has been exchanged but relative atomic dis- 
tances within each grain are still nearly unchanged. At 
t=0.02 A (fig. 6), projectile and target are still recogniz- 
able but they have traded sides and each has "expanded" : 
dissociation has begun, but only after the projectile has 
made its way through, and out of, the target. This was to 
be expected from fig. 3, which indicates, for a carbon pro- 
jectile of velocity 1000 A/ps, a range of about 50 A, much 
longer than the present target thickness (3.5 A). As to the 
target, note that its center of mass has barely moved in 
the laboratory frame up to this time. Much later, at t=0.1 
ps (fig. 7), projectile and target are shattered in 49 dis- 
tinct fragments gathered in 2 clusters of nearly equal sizes, 
and whose centers of mass are separated by about 80 A. 
Before we discuss these fragments in detail, consider first 
the energetics of the collision. Throughout the process, 
and similar to fig. 2, total kinetic and potential energy 
vary in opposite directions and total energy is conserved 
to the accuracy of the calculation. At t=0, K.E. =15816 
eV, entirely associated with the projectile. At t=0.01 ps, 
the clusters are spatially overlapping and the total gain in 



Fig. 5. Snapshot of a collision between two clusters of 
equal sizes (40 atoms) at t=0.01 ps: projectile and tar- 
get overlap. 



Fig. 6. Snapshot of a collision between two clusters of 
equal sizes (40 atoms) at t=0.02 ps: projectile now on the 
left after traversing the target, distance: 9 A, dissociation 
just begun. 
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Fig. 7. Snapshot of a collision between two clusters of 

equal sizes (40 atoms) at t=0.1 ps: nearly complete va- 
porization, 49 fragments, mean distance between clusters 
^ 80 A. Note that the scale is much larger than for fig. 4 
to 6. 

P.E. peaks at about 1400 eV, with a corresponding nega- 
tive peak in K.E. (deceleration of projectile by the poten- 
tial barrier of the target atoms). The full width of these 
peaks is 0.007 ps, which is the time spent by the projectile 
inside the target for this initial velocity. During this time, 
a fraction of the translational K.E. of the projectile is con- 
verted into internal K.E. in both clusters. This energy, in 
turn, is expended in bond dissociations and transformed 
in P.E. and K.E. of the fragments. At the end of the re- 
action, the two groups of fragments nearly equally share 
the total gain in P.E., which is ^262 eV. For 49 products 
(fig. 7), this amounts to an average of 5.3 eV per dissoci- 
ation. On the other hand, the final K.E.s of the remains 
of target and projectile are very different: the latter have 
retained ~88% of the initial K.E. while the former have 
only gained 61 eV; and the velocity of the latter's center 
of mass is still a full 940 A/ps (in the initial direction), 
while the velocity of the former's center of mass is only 60 
A/ps (nearly in the same direction). 

Similar numerical simulations were performed for 
smaller and slower projectiles. Sizes of 40 (24 C and 16 
H), 20 (12 C and 8 H), 10 (6 C and 4 H), 4 (2 C and 2 
H) atoms as well as a single C atom (Sect. 3) were used. 
The corresponding set of results is represented in fig. 8, 
where the abscissa is the initial projectile velocity and the 
ordinate is the ratio of the number of atoms in fragments 
dissocited from target and projectile, to the total number 
of atoms in both clusters: this is taken as measure of the 
dissociation efiiciency of a collision. The following can be 
noted: 

— For all particle sizes, dissociation begins at Vi~100 
A/ps (10 km/s), when each C atom of the projectiles 
carries ~6 eV of energy, just enough to break a carbon 
bond. 

— For all projectile sizes, a maximum number of disso- 
ciations is essentially achieved at 14^200 A/ps, when 
according to fig. 3, the range of a C atom into the 
target is just about equal to the thickness of target 
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Initial relative velocity (A7ps) 

Fig. 8. Dissociation efficiency of collisions: ratio of number 
of atoms in fragments dissociated from target and pro- 
jectile, to total number of atoms in both clusters; as a 
function of initial relative velocity. Crosses, plusses, up 
triangles, circles and a square are used, respectively, for 
projectiles of 40 (24 C and 16 H), 20 (12 C and 8 H), 10 
(6 C and 4 H) , 4 (2 C and 2 H) atoms as well as a single C 
atom; the target is a 40-atom cluster ( same as in fig. 4) for 
all simulations in this figure. Targets and projectiles (ex- 
cept single C atom) have same thickness, 3.5 A, parallel to 
the initial velocity and to the line joining centers of mass. 
The size of projectiles varies in the plane perpendicular to 
this direction. To help the eye, lines are drawn between 
discrete simulation results. 1000 A/ps=100km/s. 

or projectile in the present series of experiments (3.5 
A). Beyond this, unused energy is carried away by the 
remains of the projectile leaving the target. 

— The maximum number of dissociations that is achieved 
increases with the transverse size (hight) of the projec- 
tile relative to that of the target, but not much beyond 
a size ratio of ^0.5. In this series of simulations, this 
ratio corresponds to a difference of transverse sizes of 
only 10 A, or about the cross-section diameter of a car- 
bon atom suggested by the results of Sect. 3. Thus, it is 
likely that, for grain sizes much larger than those con- 
sidered here, maximum dissociation efiiciency would 
be reached only for nearly equal transverse sizes of 
projectile and target. 

— Comparison of simulation results suggests that, as far 
as dissociation efficiency is concerned, the effect of a 
projectile can be considered as the sum of the effects of 
its individual carbon atom constituents. For instance, 
for given target size and initial projectile velocity, the 
final gain in total P.E. (which is proportional to the 
number of resulting dissociations) is nealy proportional 
to the number of C atoms in the projectile. 

Drawing upon the same set of simulations, fig. 9 plots 
the number of products (= number of broken bonds) as a 
function of the final gain in total P.E. The average ratio 
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Fig. 9. The average final gain in total P.E. per dissociation 
(or broken bond) is the inverse slope of the straight line 
fit: 5.6 eV per bond. leV=23 kcal. 

of these quantities is found to be 5.6 eV per bond. On the 
other hand, the average bond energy of the model material 
may be estimated as the ratio of total initial P.E. (9227 
kcal or 401 eV) to total number of atoms (80), giving 5 eV 
per atom, which is lower than the average dissociation en- 
ergy just derived. The difference of 0.6 eV per dissociation 
indicates that the fragments also carry some vibrational 
energy, which amounts here to an average temperature 
of at most 4660 K divided by the number of atoms in a 
fragment. All this vibrational energy is but a tiny fraction 
of the initial projectile K.E.; it increases as the latter de- 
creases towards the threshold (14~100 A/ps) but always 
remains small. 

Other interesting results of these experiments involve 
the average and maximum sizes of the products of colli- 
sions (fig. 10). As the projectile initial velocity increases 
beyond threshold, it penetrates deeper into the target ( cf. 
fig. 5), dissociating more bonds and creating more frag- 
ments of decreasing size (number of atoms). This trend 
levels off, as expected, beyond the velocity which is re- 
quired for the projectile to make its way throughout the 
target, i.e. about 200 A/ps in the present instance. For 
a given initial velocity above this limit, the average and 
maximum sizes of products decrease as the transverse di- 
mension of the projectile gets closer to that of the target. 
At equality, which is the case of fig. 10, the average size of 
products is then smaller than 2 atoms, meaning that, as 
observed, there are many free C and H atoms among the 
fragments; also, the maximum size of the latter is about 10 
atoms, meaning that the fragments are molecular rather 
than granular. Within the limits of validity of the present 
treatment, there is no reason to believe that this result 
depends on the initial sizes of projectile and target. We 
are therefore led to surmise that, in most circumstances, 
vaporization is the main destruction process affecting the 
grain size distribution. The next section considers this is- 
sue in more detail. 



Fig. 10. Maximum (triangles) and average (squares) num- 
ber of atoms in fragments resulting from a collision be- 
tween clusters of equal sizes (40 atoms), as a function 
of the initial relative velocity. The fragments are essen- 
tially small molecules as soon as the range corresponding 
to the relative velocity exceeds the size of the clusters. 
1000 A/ps=100km/s. 

6. Grain destruction by collision 

From the previous sections, we can draw the following 
simplified picture. Consider a target grain of size D and 
a projectile of size d and initial velocity V towards the 
target. Let R(V) be the range of the projectile through the 
target material. Several generic cases must be considered 
for the outcome of the collision. 

1. D=d 

(a) d<R: complete vaporization of both grains (cf fig. 

7 and 8); 

(b) d>R: part of the projectile (volume v=d^R) pene- 
trates into an equal volume of the target and both 
volumes axe vaporized; formation of a crater of 
depth R and transverse size d in both grains; 

2. D>d 

(a) R<d<D: vaporized volume v= 2d^R; cratering in 
both grains; 

(b) d<R<D: v=d^(R-|-d); cratering in target, and pro- 
jectile entirely vaporized; 

(c) d<D<R: v=d^(d-|-D); complete vaporization of 
projectile and of part of target (hole of diameter 
d). 

As in fig. 8, let the destruction efficiency of the collision, 
f, be defined by the ratio of the vaporized volume v to 
the total grain volume d^+D^. The different generic cases 
can be synthesized as in fig. 11, where Vd and Vc2 are 
critical velocities respectively defined by the conditions 
R=d and R=D. Assume, for simplicity, that the range R 
is proportional to V. Then, f has simple expressions in the 
three intervals of interest (see fig. 11): 

- 10 km/s<V<Vci: f=2RdV(D3-Kd3) 
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Fig. 11. Destruction efficiency, f, of collisions between 
clusters of sizes D and d, as a function of initial relative 
velocity (cf fig. 8). Vthr~10 km/s. Vci and Vc2 are critical 
velocities respectively defined by the conditions R=d and 
R=D, where R is the penetration range relevant to the 
model material (see text); f strongly depends on relative 
size of projectile. 

— V,i<V<Vc2: f=(R+d)dV(D3+d3) 

— Ve2<V: f=d2(l+d/D)/D2(l+d3/D3). 

Hence the following cursory inferences: 

— No destruction occurs for V<10 km/s. 

— Sticking and aggregation (coalescence) of grains oc- 
cur over a small range of velocities about 10 km/s. At 
lower velocities there is quasi-elastic recoil. The case 
of a single-atom projectile stands apart because of the 
strong chemical reactivity of single atoms. 

— The collision destruction efficiency is highest for rela- 
tive velocities higher than Vc2, and only reaches 1 for 
reactants of equal sizes. 

— C-C and C-H bonds stand nearly equal chances of be- 
ing broken. 

— Whatever part of the grains is destroyed, is not frag- 
mented into intermediate-size grains, but reduced to 
molecules of a few atoms (vaporized). This seems to 
question the role of collisions in producing the famous 
grain size distribution in r""^-^ (sec also Thebault et 
al. If this is confirmed, then the origin of such 
a distribution must be sought elsewhere. 

— For V=100 km/s (or 1000 A/ps), R-80 A and com- 
plete vaporization occurs only for equal grain sizes 
smaller than 80 A. 

— Equation 2 cannot be extrapolated to higher velocities 
because the present treatment does not include elec- 
tronic interaction between projectile and target (see 
Sigmund ll969|l . However, experimental and theoretical 
results indicate ranges increasing with projectile K.E. 
to powers of 0.5 to 0.75 up to hundreds of keVs (see 
Schiott [l9f)8l Biersack and Haggmark 1980}. These re- 
sults suggest that complete vaporization of the largest 



grains in the ISM ('^ O.l/im) requires relative velocities 
of thousands of km/s. 
— For high velocities and unequal grain sizes, part of 
the big grain and all of the small grain is vaporized 
into small molecules which carry, in the form of trans- 
lational energy, most of the coUisional energy that 
was not used for dissociation. The residual grain frag- 
ment(s) do(es) not reach temperatures high enough for 
it to suffer significant structural changes as compared 
with the initial material of the reactants. This conclu- 
sion is partly dependent on the microscopic, dynamic 
approach used here, as opposed to the more usual 
macroscopic, thermodynamic treatment. The former 
implicitely assumes that all translational energy car- 
ried by fly-away fragments is lost to further dissoci- 
ation or heating, which is not the case in the latter. 
This must make a difference for the largest grains. 

In conclusion, the present atomistic, non-equilibrium, dy- 
namical approach sheds a different light on, and warrants 
some caution in, assessing the threshold for grain destruc- 
tion and the impact of collisions on grain size distribu- 
tion as well as on grain processing. This detracts noth- 
ing from the validity of the shock wave treatment in cases 
where most of the projectile energy is delivered in the bulk 
of the target by a wave travelling at a velocity of a few 
km/s. Such is the case when the projectile size and range 
in the target are much smaller than the size of the latter. 
Numerical simulation is also applicable in this situation, 
and can be very helpful in describing crater formation and 
temperature and pressure ditribution throughout the tar- 
get as a function of time (see Cleveland and Landeman 
ITTMI Aoki et al.HMZIl- 
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